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Abstract 

It is shown that the Poisson bracket with boundary terms recently 



proposed by Bering ( hep-th / 9806249 ) can be deduced from the Pois 



son bracket proposed by the present author ( |hep-th/9305133| ) if one 
omits terms free of Euler-Lagrange derivatives ("annihilation princi- 
ple"). This corresponds to another definition of the formal product 
of distributions (or, saying it in other words, to another definition of 
the pairing between 1-forms and 1-vectors in the formal variational 
calculus). We extend the formula (initially suggested by Bering only 
for the ultralocal case with constant coefficients) onto the general non- 
ultralocal brackets with coefficients depending on fields and their spa- 
tial derivatives. The lack of invariance under changes of dependent 
variables (field redefinitions) seems a drawback of this proposal. 
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1 Introduction 

Recently Bering proposed a formula for the field theory Poisson bracket 
with boundary terms which are different from those proposed earlier by the 
present author . In general the motivation for new brackets arises from the 
fact that the well-known standard field theory Poisson bracket is not applica- 
ble to many nontrivial boundary problems. In particular it does not satisfy 
the Jacobi identity The terms violating the Jacobi identity are, of course, of 
purely boundary (or divergence) nature and so can be killed by some bound- 
ary conditions which are called trivial here. The problem addressed in this 
paper as well as in publications |T], 0, |^, || consists in searching for the Pois- 
son bracket formula which exactly fulfils the Jacobi identity even before any 
boundary conditions are taken into account. 

According to our knowledge the first observation of this problem and the 
first successful attempt to solve it was made by Lewis, Marsden, Montgomery 
and Ratiu (LMMR) in their treatment of dynamics of the ideal fluid with a 
free boundary || . Both formulae proposed in |IJ and [[J are in fact only two 
different extrapolations of the same formula suggested in 0. So, it would be 
better to remind first the pioneering approach. 

In contrast to the popular view [EJ that in the Hamiltonian approach all 
functionals should be "differentiable" the variations of functionals studied by 
LMMR are not free of boundary terms 



n an 

The idea is to generalize the definition of the variational derivative by incor- 
porating the boundary contribution 
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Then LMMR in fact proposed to use for the new Poisson bracket the old 
formula 
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One immediately sees that the most dangerous term with the product of 
5-functions is absent above. In fact, to kill this term LMMR put a special 
boundary condition 

5 V F5 V G 5 V G5 V F 
5q 5p 5q 5p 

which enforces zero value for the coefficient standing before this dangerous 
product .0 Unfortunately it is not quite clear whether the Poisson bracket 
{F, G} preserves this property in general case even if the initial functionals 
F and G satisfy ©. 

Here we see the bifurcation point for the following generalizations of 
LMMR result. The idea of is that these dangerous terms with products 
of 5-functions must be omitted independently of any boundary conditions 
("annihilation principle"). Another idea advocated earlier in is to find a 
reasonable formula for these terms. 

To explain this in more detail we need first to introduce a relevant for- 
malism for treating general variations of functionals depending on arbitrary 
(but finite!) number of spatial derivatives. In the adequate mathematical 
machinery was found to be the so-called higher Eulerian operators || |9|, 10 1 . 



2 It seems interesting to mention that in recent calculations of the central charges arising 
in boundary symmetry algebras in Chern-Simons theory Q the similar term is cancelled 
automatically due to the Poisson structure itself. 
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We shall follow notations of @, (§. Einstein rule is used, i.e. we omit the 
sign of summation over repeated indices and multi-indices. 

The first variation of a general local (max| J| < oo) functional 

j f{4> A {x),4> { l\x))d n x 

h 

can be written in a form 
SF = j -^ Jj Dj5 ( p A d n x = J f A (HA)d n x = j Dj {E J A {f)5<p A ) d n x, (6) 

where in general J denotes multi-index J = (ji, ...,j n ) and 

(J) dlJl< pA , , 

4* A = ~ryf, i -K— - = Dj(j) A , \J\=ji + ...+j n . 

but in the simplest case of one- dimensional space it is just the order of spa- 
tial derivative. We have introduced also the Frechet derivative which is a 
differential operator 

f A = {LD, (7) 

The higher Eulerian operators E J A are uniquely defined by the following for- 
mula 

Here binomial coefficients for multi-indices are 

J \ fh\ fir 



KJ \kj \k 

j\ = (j\/(k\(j-k)\) H0<k<j; 
kj \ otherwise, 

and 

{-D) K = (-1)WD K . 

Both and use the full variation (Q) for the construction of the new 
Poisson brackets but in different ways. In it was proposed to start from 
the formula 

{F,G} = 5 G F = J Dj{E J A {f)5 G ^ A )d n x = J f A MA)d n x, 



n 



and to look for 8 G (f>A of such a form which fulfils the equation 

S G F = -5 F G. 

The following formula was derived in after some calculations 

{F,G} = J D J+K [E J A {f)i AB E^{g)) d n x = J Tr (j'JabQb) d n x. 
n n 
We can for the easier comparison with |l| first consider only the so-called 
ultralocal brackets then 

{4>a(x), 4>B(y)} = IabS(x, y), I A b = -Iba- 

In contrast, the proposal of [[TJ is to start with the already antisymmetric 
expression 

{F, G} = A G F — A F G — {F, G} id, 

where 

{F,G} old = J E° A (f)I AB E° B (g)d n x, 
h 

A G F = J Dj(Ei(f)A G ^ A )d n x = J f A {A G ct>A)d n x. 

n n 
Then it is possible to use the standard expression for the field variation 

Ag^a = IabE° b {9) 
and the resulting formula will be 

{F,G} = J Dj(Ei(f)I AB EUg)-Ei(g)I AB E B (f)y n x-{F,G} old . 
u 

So, it is easy to see that the last formula contains only one summation 
over multi-index J whereas formula (0) contains a double sum over J and K. 
If we omit all the terms without at least one of E° operators in this double 
sum we immediately get (Q). 

Maybe it will be of some interest to add that in ultralocal case for the 
local functionals depending on the spatial derivatives of the fields of order 
up to N, Bering's bracket involves spatial derivatives of order 3iV, whereas 
the bracket proposed in p[ involves 2N, as also the standard bracket does. 

The point of difficulty with Bering's formula seems to be the lack of in- 
variance under the changes of dependent variables (differential substitutions 
of fields). 
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2 Differential substitutions 



Let us consider the invariance properties of the field theory Poisson brackets 
under field redefinitions of the type 

<\>A -> <t>B = £bOa, Dj(f) A ), (9) 

(differential substitutions). 

If we initially have some local Poisson brackets for fields 4> A (x), i.e. 

{(j) A (x),(f) B (y)} = I AB (x)5(x,y), (10) 

where I AB — I A bDk is a differential operator of a finite order with field- 
dependent coefficients 

I% B = I A ( B (ct>c,Dj ( t> c ), (11) 
then as a result of the differential substitution (|9p we get a result 

{<f>c(x),(f)D(y)} = (€c)'a(x) (€d)' b (y)iAB(x)$(x,y)- 

To transform this expression to the form similar to fllPf ) we need a definition 
of the "adjoint" operator 

J AB (x) 5 (x,y) = J^ omt '(y)5(x,y), 

then we will have 

(£c)'a ( x ) (£d)b (y)lAB(x)6(x,y) = (£c) A (x)Iab(x) [(£d)'b] (x)5(x,y), 
and 

Icb{x) = {ic)' A {x)I A B{x) [{Cd)' b ] " adj ° mt " (*)■ 

The approach which we consider here is different from the standard one by 
the treatment of boundary (or divergence) terms. All of them should be 
preserved in the calculations. This means that we require the exact equality 

J ZaJabVB^X = J TjA^AB^Bd^X, 

n n 

without discarding any boundary (divergence) terms. In contrast the stan- 
dard approach requires only equality up to boundary terms 

ZaJabVb ~ VaJ*ab£b (mod divergences), 
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or, in notations of Appendix C, 

(Z\J\v) = (v\J*\0- 

From the last definition we get a usual relation 

Jab = (-D)k o J* a . (12) 
But the former one gives a different result: 

JUb = (-D) k o9 u JK a . (13) 



Here we use the characteristic function of the domain of integration (physical 
domain) 

1 if x G VL\ 
otherwise; 

to codify the divergences. There is an apparent relation 



edx) 



-D) K 9 a f(x)d n x = J 6 n D K f(x)d n x = J D K f(x)d n x. 



n 



So, with 8n we are able to write all the spatial integrals not as integrals over 
the physical domain but as integrals over the whole infinite coordinate space 



3 The standard bracket 

Let us first consider transformations of the standard Poisson bracket 

{F, G} old = ~J (E A (f)I AB E° B (9) ~ E° A (g)i AB E° B (f)) d n x, (14) 
n 

under field redefinitions of the type (|9|). We will use a formula 

E° A (f) = [(Zc)' A Y EUf) = 

= {Ea(Zc) - E\fa)D + E A (£c)D 2 -...) E%(J) = 

= {-l)^Ef{i c )D K E%{fl (15) 
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derived in Appendix C. Then in the integrand of (|T^) we will have the fol- 
lowing expressions 

ElU)i AB E%{g) = [(£c)'a\* E c(f)lAB [(€d)'bY E< b(9) = 

= &c(f) [[(Zc)' A Yy '1ab[{£d)' b Y E °d(9)- (16) 

But, as we discussed above, the two definitions of the "adjoint" differential 
operator are not equivalent 

mc)' A YY=(Zcy A oe u ^(Zc)' A , 

and so 

Of course, the difference has a form of divergences. Therefore we see that 
the standard field theory Poisson bracket is invariant under field redefinitions 
of the form (|j) (i.e. differential substitutions) only up to boundary terms. 
Surely, this is enough because all the other requirements (Jacobi identity, 
antisymmetry, Leibnitz rule) are also fulfilled only up to boundary terms 
and this bracket does not pretend to be adequate for nontrivial boundary 
problems. 



4 The general approach to both proposals [[Hfl 
and [@] 

Dealing with local functionals and local Poisson brackets in field theory we 
always get a general bracket of two functionals in the following form 

{F,G} = J d n x J d n yf A (x)g' B (y){Mx),My)}, 
n n 

where Frechet derivatives f' A (x), g' B (y) are differential operators (0). If we 
represent each integral not as the integral over finite domain but as an integral 
over all the infinite M™ with a characteristic function of the domain Q 

{F,G} = J 6 n (x)d n x J en(y)d n yf , A (x)g' B (y){M^,My)}, 
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then it easy to integrate by parts formally and get 

/f ftF fiC 

dn% / dn y ju ( \ A/A ( A M*), My)}, 
J d(p A (x) d<p B {y) 



1" B» 



where 

t^t = *n = (-D) K [Bu-^L ) = EXJ){-D)j9 a , 

with the the analogous formula valid for the variational derivative of G. Let 
us restrict the consideration to ultralocal case for simplicity, then 

{F,G}= [ d n x [ d n y-^I ABJ ^5(x,y), 
J J d<p A {x) o<pB{y) 



tn ran 



and in both approaches it is believed that 



{F, G}= [ d n x- 
J < 



SF , 5G 



A (x) AB 54> B (x)' 

But, of course, this expression includes products of distributions of the form 
DjO^i x D^On and here these two proposals are different: 

• in Bering's approach 

Dj9n x DkOq = 5j D k 9q + 5koF>j9q, — 5j 5ko9q, 

• in the approach of 0: 

Dj9q x D k 9q = D j+k 9q. 

Apparently it is possible to avoid these distributions completely as they serve 
only to codify divergences. This is demonstrated in publications Then 
the key transformation from the double spatial integral to the single spatial 
integral with the help of 5-function can be interpreted simply as a pairing be- 



tween 1-forms and 1-vectors of the formal variational calculus || [13], [L4| . The 



pairing defined in |] is compatible with the grading related to divergences. 
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Now if we use the above formula for Bering's pairing it is possible to derive 
the Poisson brackets in the most general (not treated in |]J) non-ultralocal 
case where Iab is a differential operator of a finite order with field dependent 
coefficients. Really, this is the same formula with hats addedQ 

{f,g} = f' A {i AB El{g))-g' A {i AB E%{f))- 

- \ (E° A (f)iABE B (g) - E A (g)I AB E B (f)) . 

Moreover, it is possible to demonstrate that the Jacobi identity is fulfilled 
for this bracket for any local operator Iab with constant coefficients . In the 
case of ultralocal Poisson brackets with the field dependent coefficients (but 
the dependence on field derivatives is excluded) in Appendix B we derive the 
following condition of the fulfilment of the Jacobi identity for the bracket 
constructed accroding to Bering proposal 

Iab,cIcd + Ida,cIcb + Ibd,cIca = 0. 

In Appendix B we also give the condition for the most general case of non- 
ultralocal brackets with coefficients depending on the spatial derivatives of 
fields also. 

5 Differential substitutions and Bering's pro- 
posal [jtj] 

Now let us consider the formula derived in the previous Section as a further 
development of the initial proposal by Bering 

{F, G} B = J f A i AB E° B (g)d n x - J g'J AB E B (f)d n x - {F, G} old , (17) 
n n 

where {F, G} \d is the standard Poisson bracket treated before. 

3 It is possible to derive a more general formula if we do not suggest that the opera- 
tor Iab should be antisymmetric with respect to the standard definition of the adjoint. 
But this formula does not fulfil the Jacobi identity even for non-ultralocal brackets with 
constant coefficients. 
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The Frechet derivative transforms under differential substitutions (Q) as 
follows (see Appendix C) 

/a = f'c (£c)a ' 

so we obtain 

fjABE B = r c (ic)' A iAB [{^)' B Y E%{g). 

It means that the first and second terms of the bracket will be invariant if 
we suppose 

Icd = {£c)a 1 ab [(€d)b\ > 

so, the old definition of the adjoint operator ( pT2|) should be used here in 
accordance with the treatment given by Bering (see Subsection (5.5) of [0]). 

Unfortunately, this bracket contains also a term {F, G} id (the standard 
Poisson bracket) with another transformation properties. As we have demon- 
strated in Section 3 this term is invariant under field redefinitions only up to 
divergences. So, taken as a whole, Bering's formula is not invariant. 



6 Differential substitutions and the bracket 
proposed in [|2j 

Let us show that in contrast the formula 



{F,G} = J Tr (j'Jabq's) d n x 



is precisely invariant under field redefinitions of the form (§). We remind that 
the trace is used here to denote the rules of composition of the differential 
operators f' A , Iab and g' B : 

Iab = IabDk, 

dg n 

Operator f' A acts on everything to the right of it, so does operator Iab, and 
operator g' B acts on everything to the left of it, i.e. acts on everything besides 
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its own coefficients 



if 



T (I'Jab'Jb) = {^jD M ^D J+K ^ M i AB D L] 

After the field redefinition we get 

{F,G} = J Tr ^ c )'J AB g' D ^ D )' B ) d n x. 
n 

So, if we use here the adjoint operator to defined by flTB| ) then it will 

act only onto g' D 

{F,G} = J Tr (f'cfoYAlAB [(Zd)'b?9'd) ^ 
n 

But according to our definitions given in Section 2 

As a result we see that this definition of the field theory Poisson bracket with 
boundary terms is exactly invariant under differential substitutions. 

In publication [I1J this invariance was demonstrated for the concrete ex- 



ample — the Ashtekar transformation 12 of the gravitational variables. 



7 Conclusion 

We considered above an interesting proposal made by Bering on the boundary 
terms in the field theory Poisson bracket. We generalize this proposal to the 
most general local Poisson brackets and find the conditions necessary to 
fulfil the Jacobi identity. According to our treatment given in more detail in 
previous publications 0, |J there are three separate ingredients of the Poisson 
bracket construction which should be revised: the differential of the local 
functional, the Poisson bivector and the pairing operation. Bering uses the 
same definition for the differential, but changes the pairing and the bivector. 
It occurs so that to change the pairing alone means to get into trouble with 
the Jacobi identity in the non-ultralocal case. 

Really, the paper ]1) suggest a lot of new ideas which deserve more dis- 
cussion. Here we only concentrated on the drawback that it seemingly had. 
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Probably, the further investigation will show whether these drawback could 
be overcome in Bering's approach. But anyhow it is absent if we use another 
formula suggested in . 

Acknowledgements 

This work was completed during the visit of the author to the International 
Center for Theoretical Physics. The author is most grateful to ICTP for 
hospitality and to Professor S. Randjbar-Daemi for the invitation. 

Appendix A. Useful relations 

Let us suppose that Iab = Iab-Dm, where I AB = I^^^), and 

I*AB = (-D)m o lg A . 

Then the following useful relations can be proved by using the technique of 
higher Eulerian operators || [HJ compiled in Lemmas 2.5 - 2.12 of 0: 

{El{g))' c {icDE* D (h)) = {-D) L 



E%(g'c{ic D E D {h))) = 

+ 
+ 

E° B (E° c (g)i CD E° D (h)) = 

+ 
+ 



d 2 



Dj(I CD E° D (h)) 



-D)l 
d 2 h 



erg 



Dj(i CD E° D (h))+ 



Dj(i* CD E° D (g))+ 



E°c(9)^D M EUh) 



■D)l 
d 2 h 



d 2 g 



Dj(i CD E° D (h))+ 



Dj(r CD E° D (g)) 



qtM 
B 



E° c (g) 
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Let also I* AB = —Iab, then from the above it follows 

E° B (g'c(ic D E D (h))) = -E% (h' c (i CD E° D (g))) 



E° B (E° c (9)ic D E D (h)) = (—D)i 



E °c(9)^D M EUh) 
d<t> B 



+ 



+ {E%{g))' c (I C nE° D (h)) - (E? B (h))' c {I CD E» D {g)). 

We will use these relations in checking the Jacobi identity in Appendix B. 

Let us illustrate our results by the less general case of ultralocal Poisson 
brackets 

Iab = Iab = —Iba, 

and suppose for the simpllicity that functions Iab depend only on fields and 
do not depend on their derivatives 



Iab = Iab{4>c)- 



Then 



and so, 



01 CD 

d4>f 



5mo8lqIcd,b, 



E° c (g)^D M EM 



I CD ,BE° c (g)E D (h). 



Appendix B. The Jacobi identity 

By using Bering's formula for the Poisson bracket we get 



{F,G} = j{f,g}d n x, 

{f,g} = f' A (iABE B (9)) - UIabEIU)) - 

- lE A (f)i AB EUg) + \E A (g)i AB E B (f), 

{{f,g},h} = {f,g} > c (ic D E D .(h))-h' c (ic D E D .({f,g}))- 

~ lE° c ({f,g})i CD E° D (h) + \El(h)i CD El({f\g}) 
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+ 
+ 
+ 
+ 
+ 



f' AC [i AB E° B (g),i CD E D (h))-g' AC (l AB E% (/), I CD Ef> D (h)) + 
!'a ((iAB)'c(icDE D (h))EU9) + iAB(E B (g)y c (i CD E° D (h))) - 
9a ((hB)'c(ic D EUh))EUf) + iAB(EUf)y c (icDE° D (h))) - 
\{El{f))' c {icDEl{h))i AB El{g) - 

lE^f)(i AB y c (i CD EUh))EUg) - 

lE M)iAB(EU9)yc(icnE° D (h)) + 
\{E\{g)y c {ic D E D {h))i AB El{f) + 
\El{g){i AB y c {ic D E D {h))E%{f) + 

lE^9)iAB(EUf)yc(ic D E D (h)) + 

h' c (l CD El{g' A {i AB El{f)) - f' A (i AB E B (g))+ 

E A (f)lA B E B (g)))) - 

\e° c (f A (iABE B (g)) - g' A (I AB E B (f)) — 



\E A {f)i AB El{g) + l -E\{g)I AB El{fyj i CD E" D {h) + 
l -E c {h)i CD El (f' A (I AB E B (g)) - g' A (I AB E B (f))- 



g cyclic permutations and applying formulae from Appendix 




Here we use notation 



f'JiB {^Ai Vb) 




{{f,g},h} + {{hj},g} + {{g,h},f} 
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-§rD L {I CD EUh)D M E%{g)) 



-9'a 



r)T M 

01 AB n (r 



D L (I CD E» D (h)D M EUf)) 



-ti 



C 1 CD 



-D) 



1 r5f M 
-^El(f)^D L (I CD E° D (h))D M E° B (g) 

1 r)/~ M 

lEO(/)^U^0,) - I^o^)^^(/) 

z 0(p c / 



ic D E° D (h) 



^E° c (h)I CD (-D) L 



\E\U)^D M El{g) - i^fo)^^ (/) 



4 L) " 



+cyclic permutation of (f,g,h). 



From the above expression it is apparent that in the case of constant coeffi- 
cients I^b the Jacobi identity is satisfied. It is straightforward to check that 
in the case of ultralocal Poisson brackets with the coefficients depending on 
the fields (but not on their spatial derivatives) we get a well-known condition 
for the fulfilment of the Jacobi identity 



Iab,cIcd + cyclic permutation of (A, B, D) = 0. 



Appendix C. The transformation rules 

Here we derive the transformation rules for Euler-Lagrange and Frechet 
derivatives under differential substitutions of fields @. 

First, let us consider the variation of an arbitrary function of the fields 

df 

94>a 
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If we use the transformed fields 



then we get 

where 
Therefore, 



df 

d<t> B 



h = (CbYa^a- 



f r A = f's° (€b)'a- 
Second, let us consider an expression 

(i|/Wa> = = E A (f)$d>A, 

where the angle brackets denote the standard integrand, defined up to diver- 
gences, and make a change of variables 



'a = £,b{4>a, Dj<p.\ 



then 



E A {f)5<p A = (l\f' s o ^ B )Wa\) = WMbMa) 

= (EUfMBy A \HA) = {8u my A y i^(/)>, 



or, 



E A (f) = my A yEUf). 

This result can be checked by more tedious but straightforward calculation 
by using formulae for Eulerian operators given in article . 
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